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Abstract: We extend the form-factors approach to the quantum Ising model at finite 
temperature. The two point function of the energy is obtained in closed form, while the 
two point function of the spin is written as a Fredholm determinant. Using the approach 
of [0, we obtain, starting directly from the continuum formulation, a set of six differen- 
tial equations satisfied by this two point function. Four of these equations involve only 
spacetime derivatives, of which three are equivalent to the equations obtained earlier in 
0,0]. In addition, we obtain two new equations involving a temperature derivative. Some 
of these results are generalized to the Ising model on the half line with a magnetic field at 
the origin. 
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1. Introduction 



The form-factors approach has been a surprisingly efficient tool for computing cor- 
relation functions in 1 + 1 quantum integrable theories at zero temperature (T) (see eg. 
01 1 HI > ) [0 j 1U ) • Recall that a quantum integrable theory in 1 + 1 is conveniently described 
in terms of particles (massive or massless), with a factorized scattering encoded in the S 
matrix (solution of the Yang-Baxter equation). Multiparticle states provide a natural basis 
for the space of states of the theory @ , @ . The form-factors are the matrix elements of the 
physical operators in this basis. In most theories, they can be obtained by solving a set 
of axioms generalizing Watson's equation Q]. Correlation functions are then expressed as 
infinite series of products of these matrix elements. While exact, re-summed expressions 
are hard to obtain, the approach is usually very satisfactory from a practical point of view: 
the series converge rather quickly allowing in particular potential comparison with 

numerical simulations. 

So far, only integrable quantum theories at T = have been considered (equivalently, 
statistical mechanics models on the plane). For the purpose of comparison with exper- 
iments however, it is crucial to obtain results at finite temperature (which corresponds 
in statistical mechanics to models on a cylinder). This seems to be a rather challenging 
subject. In [110], a related problem has been considered, that of T = but a non- vanishing 
chemical potential. The main difficulty encountered there was that the non vanishing po- 
tential changes the nature of the ground state, and to use the form-factor approach one has 
to compute matrix elements of operators between states with a finite number of particles 
per unit length. This is rather delicate; in fact, simply computing the scalar product of 
such states has been a major achievement in the study of integrable systems |TT|] , ]I| . 



In the present paper, we present a first step towards finite temperature results by 
treating in detail the quantum Ising model in one dimension (Id). Many of the difficulties 
expected in general integrable systems, disappear in that case due to the simplicity of the 
S matrix, S = — 1. Some of them remain however, especially for the spin operator which 
is non local in terms of fermions. 

The Id quantum Ising model also presents physical interest on its own. It presents 
a quantum-critical point between an ordered ground state with a broken symmetry and 
a quantum paramagnetic phase, and is the simplest system for which finite temperature 
crossovers near quantum critical points can be exactly computed. The crossovers are non- 
trivial and most of their qualitative features do generalize to higher dimensions. These 
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exact results are therefore an important testing ground for our current understanding of 
these crossovers [I2~| []T3| [H . Indeed, as was recently noted |T5[], there are a number of 



striking similarities between the crossovers of the d = 1 quantum Ising model and those 
of the d = 2 0(3) quantum rotor model (the latter is a useful model of anti-ferromagnetic 
spin fluctuations in the cuprate superconductors |Tjl,[|rjJ].) Finally, it is not unreasonable 
to hope that detailed low temperature studies of experimental realizations of the d = 1 
Ising model will emerge in the future. 

The paper is organized as follows. In sections 2 and 3 we develop a method to handle 
correlators at finite temperature. This method focuses onto excitations over the "thermal 
ground state" (rather than the usual vacuum of the theory) , and seems the most suitable to 
generalizations to truly interacting theories [IU] . The method is explained in more detail in 



Section 2, where we derive form- factors expansions for the d = 1 quantum Ising correlators 
at finite temperature. The two point function of the energy is explicitly computed, and 
a Fredholm determinant expression is obtained for the two point function of the spin. In 
Section 3, the analysis is extended to the model on the half line with a magnetic field at 
the boundary. Again, the one point function of the energy is computed, and a Fredholm 
determinant expression for the one point function of the spin is obtained. In section 4, we 
return to the bulk problem, and obtain differential equations satisfied by the two-point spin 
correlators. Earlier work had obtained differential equations for the d = 2 classical Ising 
model in the infinite plane (equivalent to the d = 1 quantum Ising model at T = 0), and 
precisely the same equations were later shown in || to also apply to the d = 1 quantum 
Ising model at non-zero T. Among our results is a set of similar equations which are 
consistent with, but slightly stronger than, those of ||, and obtained by a very different 
method. In addition, we will obtain entirely new equations which have the novel feature 
of involving derivatives with respect to temperature, and which are therefore special to 
the non-zero T case. Differential equations for correlators with temperature derivatives 
were obtained earlier in |1J for the non-relativistic dilute Bose gas, and the methodology 
of Section 4 is similar. In the conclusion, we comment on the differential equations that 
should be satisfied by the one point function of the spin with a boundary. We also put our 
results for the two-point spin correlator in the bulk in a more physical perspective. 



2. The form-factor approach to correlators of the Ising model at finite T. 

2.1. Generalities 



We describe the quantum Ising model in Id in terms of massive fermions. The asso- 
ciated creation and annihilation operators have the algebra: 

A(^)A((3 2 ) = -AUfyAOh) 
Ai(/h)Ai(lh) = -Ai(lh)A*(lh) (2.1) 
A^ X )A\^) = -A\(3 2 )A((3 1 )+2it5((3 1 -(3 2 ). 

In these expressions, the /Vs are the usual rapidity variables parameterizing momentum 
and energy as p(/3) = m sinh /?, e(/3) = m cosh (3. 

Let us now consider the question of evaluating the two point function of some operator 
O in the Ising model at finite temperature T. The space coordinate is x, — oo < x < oo, 
and we call t the time coordinate. Formally one has: 

< O(«, t )O(0,0) >= i E *- B * IT < M ° {X 'T, >K t',!T m> - (2-2) 

Z f-* < ib lib >< ib 'lib '> 

where the states \ib > are multiparticle states, the matrix elements of O are the so- 
called form-factors, and Z is the partition function which formally corresponds to Z = 

There are divergences in the expression (^[^) that are a result of working with multi- 
particle states that are defined in infinite volume (the x-direction) . These infinite volume 
divergences manifest themselves in the norms of the states, e.g. < /3\/3 >= 2n5(0). These 
can be regulated by letting the cylinder have a finite length L and letting L —>■ oo. The 
5(0) are then regularized to cosh (3. With this explicit regularization, it is natural to 
consider a thermodynamic approach to the expression ( |2.2| ), which we describe in the re- 
mainder of this section. In this free field situation, it is possible to deal with the expression 
( p.2|) directly, and this provides a check on the thermodynamic approach; this is presented 
in Appendix A. 

2.2. Thermodynamic approach 

If we regularize the norms as explained above, the sum over intermediate states is a 
discrete sum over states with allowed momenta. On general grounds, we expect the parti- 
tion function to go like exp(Lf) where / is the dimensionless free energy per unit length. 
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That Z behaves exponentially means that the sum is dominated by terms where the energy 
goes like L, the size of the system. In other words, the sum is dominated by multiparticle 
states that have a number of particles proportional to the size of the system L. These 
states we call macroscopic, and for our purpose they will be completely characterized by 
non- vanishing densities of particles (per unit length) as L — > oo. This situation is to be 
contrasted with computations at T = 0, where one considers only microscopic states, that 
have a finite number of particles as L — > oo. 

Since microscopic states give a negligible contribution to ( |2.2|) for T ^ 0, we can 
replace by a functional integral : 



[dp], (2-3) 



where p denotes the density of particles per unit length (the exact measure to put in 
deserves some discussion, but it will not matter in what follows). The denominator and 
numerator become respectively 



[dp]e 



-0{E([p])+S([p\)} 



J y < w > 

Here E[p] = mL J d9p(9) cosh#, the entropy follows from Stirling's formula S = L J dd[(p+ 
p h ) ln(p + p h ) — p In p — p h In p h ] , where p h is the density of holes, p the density of particles, 
with p((3) + p h ((3) = cosh (3. The point now is that the argument of the numerator 
and the denominator ( p.4|) differ only by the correlator, which is of order one. In the 
limit L — >• oo, both are dominated by the same saddle point, characterized by p/p h = 
exp(— m cosh/9/T). From now on we reserve the notation p for the saddle point density. 
Then the exponential contributions of the numerator and denominator cancel out, and we 
get the simple result (a similar result appears in a different context in 0]) : 

^nt *wnn» 1 \- < 0r|O(s, W >< ^|0(O,O)|Or > (t) 
< t)O{0, >= < Qt|Qt > ^ , (2.5) 

where we denote by \0t > any multiparticle state that is characterized by the macroscopic 
density [p] (which one does not matter). We will refer to this state, a bit incorrectly, as 
the "thermal ground state" . 

Now, the remaining task is to evaluate the matrix elements. There are several ways 
to proceed, - here, we present what seems the most physical approach, and probably the 
most suitable to generalizations. 
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2.3. Excitations over the thermal ground state 

The idea is not to use the usual form- factors, which correspond to excitations above 
the vacuum |0 >, but rather to introduce new form- factors appropriate for the excitations 
above \0t >• Consider therefore a multiparticle state that realizes the density p. The 
excitations above this state are of two kinds: one can add a particle, or one can create 
a hole. We will take into account these two possibilities by introducing two types of 
"thermal" particles, characterized by a label e = ±1. Observe that the sets of rapidities 
where a thermal particle with e = 1 and e = — 1 can be created are separate. The 
creation/annihilation operators of thermal particles with the same e satisfy ( |2.1|) while 
operators with different e anticommute. The excitations therefore still have factorized 
scattering with S = — 1. We now define the form- factors of an operator O as 

F((3l, P 2 , — ,/3 n )ei,...,e„ = , = £ "'-' £l <(3n...Pl\O\0 T >, (2.6) 

where the subscript T is implicit for multiparticle states. These form-factors have to satisfy 
relations that are very similar to the ones for excitations over the ground state. One has 

F(Pl, @2, Pn)ei,...,e n = ~F(P 2 , Pi, Pn)e 2 ,e 1 ,...,e n , (2.7) 

from the fact that S = — 1. CPT invariance gives : 

F(l3 1 ,l32,...,Pn + 2m) ei ,..., £n =F(0 n ,p 1 ,...,0 n - 1 ) £n , £li ... i£n _ 1 . (2.8) 

The third axiom is slightly more complicated. First, one expects as usual an annihilation 
pole at Pi = Pj — Z7T when q = €j. In addition, remembering the origin of \0t >, and the 
fact that a hole and a particle can annihilate each other, there should also be a pole at 
Pi = Pj when 6i = —€j. This leads to 

Res^ 2 _ in F( n \p 1 ,p 2 ,...,p n ) ei>ei> ... en =F^- 2 \p 3 ,...,p n ) £3 _ £n 

(2.9) 

Res f3l=f32 F^(p l ,p 2 ,...,p n ) £u _ £l _ £n =F^- 2 )(/3 3 ,...,/3 n ) £3 ,..., e „ 

There are two minimal solutions to these equations. The one corresponding to the 
energy operator is 

F(p 1 ,p 2 ) ++ =msmh(^-^) 

F(p 1 ,p 2 ) + - = im cosh( — ^ 2 ) 

P1-P2 (2 ' 10) 
F(Pi,p 2 ) 1- = — imcosh(- 



F{p x ,p 2 )— =msinh( 
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P1-P2 



The solution corresponding to the spin a and disorder operators fi is 

m, /?nk,..,e„ = i n/2 n tanh i^r) ^ > ( 2 - n ) 

where n is even for a, with < OtMOt > / < 0t|0t >= 1, and n is odd for fx. Note that 
we have chosen to work in the ordered phase here, and that the form-factors approach 
requires us to consider the subtracted energy operator (of vanishing one-point function), 
while the spin is not subtracted. In (|2.10| ) as well as ( |2.11|) , the overall normalization is 



largely arbitrary. 

In this free field situation, the particle and hole form factors are deduced from the 
usual ones for particles only by shifting the hole rapidities by iir. This hints at a more 
direct way of obtaining expressions for the correlators, which is explained in appendix A. 

Knowing the form-factors, the two point function of the operator O can finally be 
computed. A little subtlety appears in the sum over intermediate states. While multi- 
thermal particle states are normalized as usual, their allowed rapidities are distributed 
according to the filling fractions 

■^(^) = l + g-emcosh/3/T- ( 2 - 12 ) 

One finds therefore 

< O(x,t)O(0,0) >=f^^Yl I f[!.^f e M)e~ l£%[mtCOShPz ~ mxS ' inh0l] } 



n=0 £i >/ ~ txJ i=l 

|2 



(2.13) 



x |F(/? l5 ...,/3, 



• j Mnjex,...,e n I 



It is not difficult to verify that the above expression satisfies the Martin-Schwinger- 
Kubo relation: 

< 0(x, t - i/T)O(0) >=< O(0)O(x, t) > . (2.14) 

2.4- The energy correlations 

As a check, it is useful to consider the conformal limit of the energy two-point function, 
which can be compared with standard results. We obtain the conformal correlator by the 
limit 

P = ±P ±d, Po^oo 

m a (2.15) 
m — > 0, — e Po — > fj,. 
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Here (jl is a parameter adjusting the normalization. The spectrum separates into left and 
right moving particles with e = fie e = ±p. In this limit, only the left-right and right-left 
form factors will be non zero. They are simply : 

F(0i,fc)J£=vie<*+*»Va 

F(e 1 ,e 2 )*l = -i^ +e ^/ 2 (2 ' 16) 

The two point function of the energy factorizes as the product of a right and a left moving 
part. For the right moving part we find : 



-oo ^ 



e 2z ^ e ]e e 



oc 
oo 



du 



e —2uz R ^2uz R 



+ 



1 + e~ u / T 1 + e u l 



(2.17) 



2 sin(27rTzi?) ' 

where we have set zr = (t — ix)/2 (and similarly zl = (t + ix)/2). This is the well 
known conformal result for the fermions correlator, and multiplying by its conjugate gives 
of course the correct energy correlator. 

In the massive case, the holomorphic and anti-holomorphic parts will not be decoupled 
and the solution will be more involved. Equation ( |2.13| ) gives, restricting to x = for 
simplicity : 

^— f- < e 0,t e 0,0 >= 

> e imt(u 1 + l/u 1 )/2 e imt(u 2 + l/u 2 )/2 / y | \ 2 

„ rfu l rfu 2 " - em(ui + i /ui)/2 T 1 + e m( M2 + l/« 2 )/2T ~ ~J 

/ y \ 2 e »mt(ui+l/ui)/2 e -imt(u 2 +l/u 2 )/2 / -|_ ^ \ 2 

+ W~^ 2 / + 1 + e m(«i + l/«i)/2T 1 + e -m( M2 + l/« 2 )/2T + —) ( 2 - 18 ) 

g— imt(ui + l/ui)/2 e imt(u2 + l/u 2 ) /2 / y 1 N 2 

+ 1 + e -m(ui + l/ui)/2T 1 + e m( ll2 + l/« 2 )/2T + ~ 

e —imt(u 1 + l/u 1 )/2 e —imt(u 1 + l/u2)/2 / y y 

+ ' 



I _|_ e -m(ui+l/ui)/2T ]_ _|_ e -m(wi+l/wi)/2T \ Ul M2 



where we have introduced the variable u = e@ . The integrated result is given by : 

(2n) 2 , s , . 

< e(0,i)e(0,0) >= 



J^(-l) n [Kx(imb + nm/T) + K x (nm/T + m/T - imt)} ? ^ 
- I J^(-l) n [K (imt + nm/T) - K (nm/T + m/T - imt)} 



,n=0 



Here K are Bessel functions of the second kind. It reproduces the conformal result in the 
limit m —>■ 0. In the limit T — > it also reproduces the result of [16[. 



2.5. The spin- disorder correlations 

The conformal spin correlator cannot be naively recovered by the limit ( |2.15| ), the 
integrals being infrared divergent. In fact, the understanding of the spin correlator is 
considerably more involved, and will be the subject of most of this paper. 

Let us continue to Euclidean space space t — ► —it, and define z = (t — ix)/2, z = 
(t + ix)/2. Define 

r± = (a(z,z)a(0)} ± (fi(z,z)fi(0)). (2.20) 
It is convenient to express everything in terms of the variable u: 



u = e 



(2.21) 



Then, the formula ( |2.13| ) leads to 



E 



:±i 



E 



du\ duN 



N 



n^( 



Ui e 



-£i(mzUi-\-mz/ui) 



n 

i<j 



Ui — u, 



Ui + Un 



1 = 1 



(2.22) 



where it should be clear from the context that the thermal filling factors have changed 
their definition from ( 2.12 ) to one appropriate for an energy expressed in terms of the u 
variable 



/«(«) 



1 



1 _|_ e -ero(u+l/u)/2T ' 

At T = ( |2.22|) reduces to well known expressions (see e.g. \\17\\ , |T8| ). 



(2.23) 



The poles in (|2.22| ), which are due to particle hole annihilation, are now clearly seen 
as an artifact of the continuum limit. Consider the case where L is finite and the particles 
have quantized momenta. Then, a particle and a hole cannot have the same rapidity: a 
state is either filled or empty. When replacing discrete sums by integrals, this exclusion 
disappeared, and has to be reinstated by hand. The prescription will be to move slightly 
one of the rapidities off the real axis. This is what will be understood in the following. 

The functions r± can be expressed as a Fredholm determinant of a 2 x 2 matrix of 
integral operators 

r ± = Det(l±W), (2.24) 



(2.25) 



where 

W(uv)-( W++M W +-^ v ^ 

The two dimensional vector space structure arises from the e indices distinguishing particles 
and holes. Extending the Fredholm theory to this case, one has by definition: 

^ (±1)* 

N=0 ei"-ejv=± 



Detfli W) 



/"OO 

/ dui ■ • ■ dujy det{W eu 



{Ui,Uj)}, 



[2.26) 



where det{W £ii<E . (ui, Uj)} is an ordinary determinant of the finite N x N matrix with 
W ei>e . (ui, Uj) as entries: 



det {W euej (ui, uj)} 



W ClCl (ui,ui) W eie2 (u 1 ,u 2 ) 

Wc a ci( u 2,Ml) 



W 6N£l (u N ,ui] 



W, 



ei ,ejv 



[u 1 , U N ) 



W FoF „(u 2 , u N ) 



[2.27) 



Let { - } for 1 < i,j < n denote the n x n matrix with l/(ui + Vj) as entries. Then 
using the identity 



det{ 



1 , _Ui <j (ui-u j )(v i 



Ui + V 



one finds that the sum ( |2.22| ) can be reproduced with the following choice of kernel: 



(2.28) 



where 



W(u,v) 

e+(u) = 
e_(w) = 



e_|_ (u)e_|_ (v) e_|_ (u)e_ (v) 

u-\-v u—v 

e_(u)e+(u) e—(u)e—(v) J ' 

u—v u-\-v 



(2.29) 



f+(u) ( m 

exp — — (zu + z/u) 

n V I 

/_(«) / TO 

exp y + — (zu + z/u) 



(2.30) 



7T 
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In principle, all of our subsequent results for r± can be derived from the above for- 
mulation in terms of Fredholm determinants. However, the 2x2 matrix structure makes 
the manipulations quite cumbersome, and the following trick for dispensing with the 2x2 
structure is quite useful. Notice that in (|2.22| ) the double sum over particles and holes can 
be traded for an integral over u running over the whole real axis 



N=0 



(±1 



M 



30 dm 

oo 27r l«l 

n 



du 



N 



2tt\un\ 



N 



Ui e 



-(mzui-\-mz/ui) 



i<j 



Uj 



u 



Ui + u 



(2.31) 



Here, we introduced 



Now, one has 



where K is a scalar kernel 



with 



/(«) 



1 



1 _|_ e -m(u+l/u)/2T ' 

r± = Det(l ± K), 

e{u)e{y) 



e(u 



K(u, v) 



/(u)sgn(u) 



u + v 



( m r -is 

exp y~~^\ zu + z / u ) 



(2.32) 
(2.33) 
(2.34) 

(2.35) 



7T v 2 

and it is implied that the Fredholm determinant in ( |2.33| ) is taken over functions in the 
domain — oo < u < oo. There is a non-analyticity in the kernel at u = 0, and any subtleties 
associated with this have to be resolved by regarding the above formalism simply as a 
convenient shorthand for the better-defined 2x2 matrix integral operators in ( |2.24j ). In 
practice, we will find that one can simply ignore the singular behavior at u = 0: particles (or 
holes) with small u, have a very large energy (~ 1/w) and their contribution is exponentially 
suppressed. 

In Section 4 we will derive some differential equations that are obeyed by ( |2.31| ). 



3. Ising model with a boundary and a temperature. 

It is also interesting to include a finite temperature in the quantum Ising model on 
the half line x G [— oo, 0] with a boundary magnetic field. This provides a good example 
of the convenience of the thermal approach. 
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It will be useful to think of this problem from an Euclidean point of view, where the 
theory at finite temperature is interpreted as a theory on a cylinder of radius R = 1/T. 
One can then consider both possible directions as imaginary time direction; we will refer to 
the picture where the hamiltonian evolves in the x direction along the length of the cylinder 
as the L channel, and the crossed channel where the evolution is along the circumference 
of the cylinder as the R-channel. 

At T = 0, the boundary problem is very conveniently addressed in the L-channel. In 
this picture the boundary interactions are contained in a boundary state \B > at x=0. 
Introduce the quantity |19| 

K + cosh (3 



R(P) = itanh(/3/2)- 



K 



cosh (3 



(3.1) 



which is related to the reflection matrix by R(j3) = R(^ — 0). The constant k is related 
to the boundary magnetic field h by k = 1 — h 2 /2m. One has then [l!| 

1 d(3 



\B >= exp 



4n 



R{(3)A\-(3)A\(3) 



> 



(3.2) 



Suppose now that we include a temperature. As we have seen in section 2, in the 
crossed R-channel, the effect of this temperature can be absorbed by defining a dressed 
theory, with two sorts of particles, new form-factors, together with an integration metric 
/ e . If one formally carries out the same arguments as in |19f| , one ends up with a new 
boundary state that is very similar to (|3.2|), for the dressed theory (that is in particular 
the theory involving dressed form- factors) : 



\Bt > = exp 



(4 



P-ir,) Al(-(3)A\(P) 



|0t>, 



(3.3) 



where 



R+((3) = R((3), R-(f3) = R(/3 - in), (3.4) 

and the A\ operators as the same as in section 2. In ( |3.3|) , r\ > is a shift in the imaginary 
rapidity direction necessary to suppress the poles, and ensure the right T — > limit. 

For completeness, formula (|3.3| ) is proven in appendix B using a direct approach. 

As an application, let us look at the energy in the Ising model with boundary and 
temperature. Previous formula give 



4n 



(3.5) 



e 2mxecosh/3/ e /3 _ e ~(3\ 
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Observe that the integrals involving holes with e = —1 in (|3]5|) converge only when 
2\x\ < rj/T, which implies T\x\ « 1. The correlators are better defined by changing 
the contour of integration, which is actually equivalent to a computation in the crossed 
channel. For brevity, we shall not do so here. 

We can similarly consider the spin correlator. One finds 



< OtW\ - Pi, Pi, ■ ■ ; ~Pn, Pn >e 1 ,...,e n = 

n 

H) n JJtanhA JJtanh 

i=l i<j 

From this and ( |3.3|) it follows that 

< T \a(x)\B T > = 



7 ' tanh 



Pi-P.^"^ 



Pi + Pi^ 2eies 



(3.6) 



n=0 



jj - it] e 



3 2mie; cosh (3i 



(3.7) 



x 



| J tanh Pi Y\ tanh 



i=l 



i<j 



1, l;^ 2 " J 



tanh 



Pi + Pj 



2e i e j 



The reality of the one point function of a can be checked on ([T^]) by observing that 
R* e (P) = R € (-p),R* e {P) = R e (-P). 

Once again it is convenient to re-express ( |3.7|) into a simpler form by introducing the 
variable u = and trading the sum over e for an extended integral of u over the whole 
real axis. One finds then 



< OtW(x)\B t >= 



n=0 



TV. 



n 



n 



i<j 



Ui — Uj UiUj — 1 
Ui + Uj UiUj + 1 



-i 2 



(3.J 



where f(u) is as before (|2.23j ). 

As was done for zero temperature in | 20[| , one can express (|3.8| ) as a Fredholm deter- 
minant. One has the identity, 



n 

Kj 



Ui — Uj UiUj — 1 
Ui + Uj UiUj + 1 



n 



Hi - Uj 



i<j xr j 



= det 



(3.9) 
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where fx = (u + l/u)/2. Thus, 

< T \a(x)\B T > = Det(l + V) 



oo poo / \ 

^ / I ~\\dui J detV(ui,Uj), 



where V is the scalar kernel: 

e(u)e(u') 



V(u,u') 



u + l/u + u' + l/u' 



o(nl \_ i 1 ( U - 1 / u )(- iR ( u ))f( i / u ) mx(u+l/u)/2 

eW ~\i 2n\u\ 



(3.10) 



(3.11) 



Using the techniques of this paper it is also possible to compute correlation functions 
in boundary field theory where the finite size is in the ^-direction, i.e. on a infinitely long, 
finite width strip. This is described in Appendix C. 



4. Differential Equations for Bulk Correlators 

In this section we return to a study of the bulk two-point functions of spin and disorder 
operators at finite temperature. Recall that in Section 2 we had obtained an expression for 
these as a Fredholm determinant of a scalar kernel K defined in ( f2.34j ) . In this section we 
will use this expression to derive a set of non-linear, integrable, partial differential equations 
for the two-point functions. The general method we follow is that reviewed in Q for the 
case of the non-relativistic Bose gas, although we shall take a somewhat more algebraic 
point of view. The computations are somewhat involved, and readers not interested in the 
methodology can skip ahead to Section 4.5, where we give a self-contained presentation of 
the main results. 

4-1. The Resolvent 



The first step in deriving differential equations for ( |2.34| ) is to find the "resolvent" for 



the kernel K. We will need the resolvents H± defined by 

(1-R±)(1±K) = 1. (4.1) 
The resolvents R± are also integral operators, e.g. 

/"OO 

(RK){u,v)= I dwR(u,w)K(w,v). (4.2) 
13 



The resolvent for a kernel of the type ( |2.34|) was described in [21]. Let [V±] denote the 
class of scalar integral operators of the type e(u)e' (v) / (u ± v). As explained in[^TJ, these 
operators form a group with a Z2 graded multiplication law: 



[V ( _ } ] x = [V (+) ] x [V- ( _)] = [V {+) ], [V {+) ] x [V (+) ] = [V(_)]. (4.3) 

It turns out that R± has both [V+], [V_] pieces, so let us write 

R± = H±F, Ke[V-],Fe\V + ], (4.4) 

which turns out to be consistent with the Z<i gradings. Define functions g±{u) as follows: 

(1-R±)e = <7±. (4.5) 

Above, (Re)(u) is shorthand for J dvH(u,v)e(v). Then the definition ( f4.1|) and the fact 
that K T = K, R T = R also give us the relations 

(1 ± K)g± 
e(l-R±) 
<7±(1±K) 



e 

9± 
e 



(4.6) 



Then the resolvent can be expressed in terms of the g functions: 

H ( M ' v ^ = of, _ ,A (9-(u)g+(v) - g+(u)g-(v)) 

Zi\ Lb U I , 

1 ! ' (4-7) 
F ( M ' u ) = of 1 \ (9-( u )9+(v) + g+{u)g-(v)) . 

It is evident that the functions e(u) in K have the same z, z dependence as for the zero 
temperature case considered in ||21]| , thus the sinh-Gordon differential equations involving 



only z,z derivatives obtained there continue to hold at finite temperature. See also ||22|| . 
For completeness, these results will be included in the sequel. 

4-2. Intertwining Relations 

A central property of the integral operators K and R± is a set of intertwining rela- 
tions between them and a family of diagonal operators. Essentially all of the differential 
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equations and constraints that we shall derive in subsequent sections are a simple con- 
sequence of these intertwining relations. These relations constitute a succinct algebraic 
encapsulation of the integrability of the K kernel. 

The simplest of the intertwining relations are with the diagonal operators d n (u, v) (n 
integer) which multiply the operand by u n when u = v; explicitly the d n are defined by 

/oo 
dw 1 dw 2 Oi(u, wi)d n (wi, w 2 )02(w 2l v) 

(4.8) 

dwO\(u, w)w n 02(w, v) 



oo 



where Oi and O2 are arbitrary operators. Now, using the definition (|2.29| ) , we can easily 
obtain our first set of intertwining relations: 

(diK + Kdi)(u, v) = e(u)e(v) 

/ - VJ w , e(u)e(v) ( 4 - 9 ) 

(d_iK + Kd_! )(u, v) = -^^^ 

U V 

Notice that the right-hand-sides are all projection operators i.e. they are all products of 
functions of u and v. The intertwining of d n for general n with the K-t now follows by 
repeated application of the above and the identity d n + m = d n d m . One can also check 
that the operator did_i = 1 commutes with the K_|-. 

In a similar manner we can obtain the intertwining of R± with the d n by using fl4.4j ) 
and (Q: 

(diR± - R T di)(u,v) = ±g^(u)g±(v) 

(d_ lR± - H^d-ocu.v) = ± ^M^±M (4 ' 10) 

U V 

Again, the right-hand-sides are projectors. 

We now introduce our second set of diagonal operators: d z and d^, which take deriva- 
tives with respect to z and z respectively. They are defined by 



/oo 
dw0 1 (u,w)d z 2 (w,v) 
-00 



-00 



/oo 
dwOi(u, w)d z -02(w 1 v) 
-00 



(4.11) 



Notice that the derivatives act on everything to their right i.e. on the O2 and on all 
operands the left- hand-side may act on. Again, using the defining relations ( |2.29| ), ( |2.35| ), 
we can easily obtain the intertwining of the d z , d^ with K: 

(d 2 K - Kd z )(u,v) = -—e(u)e(v) 

( \ ( \ ( 4 - 12 ) 
(d-K - Kd-) (u, v) = - — 

2 u v 
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As usual, the right-hand-sides are projectors. The intertwining of the H± with the d z , d^ 
can be deduced from (|4.12| ) and (|4.1| ), but as we shall not need it in this paper, we will 
refrain from displaying it. 

Finally, we will need the intertwining of K± with an operator involving temperature 
derivatives. We shall follow the idea of Jl|, and choose an operator which commutes with 
the thermal factor ( |2.23| ). The appropriate operator, is defined by 

(Oid T 2 )(«,v) = J dw0 1 {u,w)(^-^T^+{w^^w^0 2 {w,v). 

(4.13) 

Again all derivatives act on all operators to their right. This particular combination of 
derivatives was chosen because 

u--) T-^— +( u + -) u-^-) O ((« + l/u)/T) = (4.14) 



u J dT \ u J du 

where O(y) is an arbitrary function of the combination y = (u + l/u)/T alone. This is 
precisely the combination in which T occurs in K±, and (|4.14| ) raises the possibility that 



the intertwining of with K-t- may be simple. It turns out that the result is somewhat 
more complicated than those obtained earlier, but is nevertheless a sum of a small number 
of projection operators. A lengthy, but ultimately straightforward, computation gave us 

(d T K + Kd r )(«, v) = e(u)T^-e(v) - e M T JL e M _ e{u)e{v) 

° U , V _ \ (4.15) 



-{u-v) U + -y-2 e(u)e(v) 

All derivatives, including the ones appearing on the right-hand-sides, act on all operators to 
their right. This equality holds in the sense that if both sides are multiplied by an arbitrary 
function of v,T,z,z, and v is then integrated over, the results are equal up to a surface 
term which is presumed to vanish. Notice especially the location of the T derivatives on 
the right-hand-side - this leads to the most compact form of the above equation, and is 
also the most convenient in subsequent manipulations; when the above equation acts on 
the g function, this form immediately gives T derivatives of the "potentials" which we will 
introduce in the following subsection. 

The identities ( |4.9| ), ( ^.10| ), Q4.12| ), and ( |4. 15[ ) are the main results of this section. Our 



results in the remainder of the paper for the relations obeyed by r± (except for those in 
Appendix D) follow by simple, repeated applications of these intertwining relations. 
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4-3. Definition of the potentials 

We need to introduce some additional formalism and notation before turning to the 
study of the differential equations. 

We consider here the so-called "potentials" , which are the set of possible scalar overlaps 
among the e and g functions introduced in Section 4.1. These occur, for instance, when 



we consider derivatives of the r± correlators. Thus from (|2.33|) 

d 2 lnr + = d z lnDet(l + K) 

= Tr((l + K)- 1 c> 2 K) 



m- ,, , „ . . (4.16) 
-Tr((e(l-R + )e) 1 ; 

m 

•— Tr(e# + ). 



In the last step we used the defining relation ( [4.5| ). The expression Tr (eg+) is one of the 
potentials; notice that the variable u has been integrated over, and the potential is an 
implicit function of z, z, and T alone. We will also need potentials defined by inserting 
various powers of the variable u in the scalar product: 

/oo 
duu n e{u)g+{u). (4.17) 
-oo 

To make things compact, we introduce some notation here for the potentials that turn up 
in our analysis. Our main results will be expressed in terms of the following potentials 

b± =Tr(e<7±) 

c ± = Tv(ed_ ig± ) (4.18) 
d± = (l±c T )Tr(ed_ 2 # ± ) 

The reason for the additional factors in the definitions of the d± potentials will become 
clear shortly. Not all of these potentials are independent; the intertwining relations 
and ( |4 . 1 0| ) leads to a constraint between them: 



c+ = Tr (ed_i0+) = Tr (ed_i(l - R+)e) from(|45D 

= Tr (e(l - R_)d_ie) - c+c_ from Q4T0|) (4.19) 
= Tr (<7_d_ie) — c_|_c_ from([4.6|) , 

which using (|4.17|) becomes finally 

c_-c + =c + c-. (4.20) 
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We solve this constraint by henceforth parameterizing 



c+ = 1 - 



c_ = e 



(4.21) 



Let us also generalize (|4.16|) , and tabulate all the derivatives of the r±: 



d z lnr± = T^b± 
<9jlnr ± = T^d T 



(4.22) 



Our choices in the definitions of the d± earlier were made so that the above equations 
came out in this symmetrical form. 

In our subsequent analysis, some additional potentials also arise at intermediate stages. 
To handle these, it is convenient to introduce here some additional notation 



Both ai, a 2 will drop out of our final results. 

4-4- Differential Equations 

We will first consider the differential equations which involve only derivatives with 
respect to z and z. These will turn out to be closely related to those derived earlier by 
Perk et. al. by a very different computation on a lattice model. Then we will turn to a 
new set of equations involving temperature derivatives. 

4-4-1 Spacetime derivatives 

We begin by turning the intertwining relations ( |4.12|) into differential equations for 
the / functions. Begin with the equation (1 + ~K)g + = e ( (|4.6| )) and act on both sides with 
the operator (1 — R + )d z . Simplifying this and related equations by repeated use of all the 
intertwining relations gives us 



These are a set of equations in the Lax form, and integrability condition d z d- = d-d z leads 
to differential equations among the potentials alone. However, a richer set of equations is 



ai = Tr(ed 1 ( 7+ ) + Tr(ed 1 (7_) 

a 2 = e^Tr (ed_ 3 #+) + e^Tr (ed_ 3 #_) . 



(4.23) 




(4.24) 
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obtained by acting on the definitions of the potentials ( fll8|) by d 2 and d 2 , using ( |4.24j) , 
and the derivatives of the e functions which follow from their definition in ( |2.35| ) . This 
gives us 

77? 

0,0 =y (&+ + &_) 
777 

dzd>= —(d++d-) 

Y 2 v + ; (4.25) 

c\b± = me^ sinh</> 
d z d± = roe T ^ sinh</>. 

With these equations in hand, we can now express ([4.24 ) in a more symmetrical form. 
First define the functions h±(u) by 

h±( u ) =e^ /2 g ± (u). (4.26) 

Note that e(u), g±(u), and h±(u) are the only functions which depend upon the spectral 
parameter u; all other functions are independent of u. Now it is easy to show from ( f4.24| ) 
and ( |E25D that 



(h^ 
\h-J 


)(«) 


-\{ 


\h-j 


)(«) 





—d z (f) mue ^\ / h + (u 
ue* d z (b ) 1 h-(u, , 

\ , \ / / x ( 4 - 27 ) 
9z0 me^/u\ I h + (u) 



me v lu —c\4> J \h-(u 



From (|4.22|) and ( |4.25|) , it can be shown that the second derivatives of the r functions 



are characterized by the compact equations: 

d z d z -lnr± = =f — e ±( ^ sinh</> 

J ( 4 - 28 ) 

d z c\(f) = — sinh2^>. 
2 

The second of these equations is the compatibility condition d 2 d- = d 2 d 2 of ( |4.27| ). 

Another simple consequence of ( 4.25| ) follows from combining it with the definitions 



( |4.22| ) ; one notices easily that 

d z (Hr-/r+) -<P)=0 0* (ln(r_/r + ) - </>) = 0, (4.29) 

which implies that </> = ln(r_/r_|_)+ a dimensionless constant dependent only on the ratio 
T/m. In Appendix D we show that the integration constant is in fact also independent of 
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T/m, and then show by an explicit computation at T = that its value is zero. Therefore 
we have 

e 4, = — . (4.30) 
r + 

It is also useful to consider differential equations satisfied by the a\^, potentials defined 
in ( f4.23|) ; as we will now show, these equations lead to additional differential equations 
involving only the b± and d±. Using the same methods as those used for fl4.25| ) we can 
obtain 

Tfi 

d z b± = -— (ai =f b±(b+ + &_)) 

Tfi 

dzd± = (a 2 ± d±(d+ + d-)) 

2 V ±V + " (4.31) 

dgai = —m(b + e~^ + 6_e^) cosh</> 
d z a 2 = —m(d + e 4 + d-e~ 4 ) cosh0 

The last two equations above are in fact not new; we can solve the first two for a\ 2 and 
then they follow from applications of ( |4.25| ). We can also eliminate between the first 
two equations of ( J4.31| ) and obtain 



d z b+-d z b- = ^(b+ + b-) 2 

m 2 ( 4 - 32 ) 

d z d + - djd- = -— (d + + d_) 

We will not use the a\2 any more in this paper, although we will implicitly use the first 
two equations in ( |4.31| ) in the derivation of some equations in the next section. 

4-4-2 Temperature derivatives 

We will now follow the same procedure as in 4.4.1, but with the operator replacing 
the d z and operators. 

First, it is useful to notice that 

i). (4.33) 

We can now get the derivatives of the / functions by acting on the (1 ± K)#± = e 
equations by (1 — R T )dT on the left; this gives 

yfew) <4 - 34) 
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m 












~2 







where _ 

. T _ . muz , . mz . , , . 

A = — dT( p + + 6_) + — (d_ + rf+ 

w 2 2u 

S+ = (1 + - Td T ) b+ -^(u 2 + l) + ^ (u 2 + e 2 *) (4.35) 

B- = - (1 + ^ - Td T ) 6_ - ^ (u 2 + 1) + ^| (u 2 + e" 2 ^) ; 

we have used ( f4.31| ) to eliminate the a± potential which arises in the above derivation. We 
will write this equation in a symmetrical form momentarily. 

First, we derive equations relating the T derivatives of the potentials; we will do this 
by evaluating 

r &(d_ 1 d T (e(7 ± )) (4.36) 

in two different ways. For the first we use dx(eg±) = g±dxe + edxg± and evaluate the 
right hand side using ( |4.33[ ) and ( [4.34]) ; alternatively, we can write, using the definition 

Tr (d_id T (e0±)) = (Td T - 1) Tr (e(l - d_ 2 )g ± ) , (4.37) 

where we have integrated by parts over the u integral. Comparing the results of these two 
evaluations, we obtain 

(1 + zd z - Td T ) b± - e ±(j} (l + zd*- Td T ) d± = m(z - z) sinh <j> (4.38) 

Again, we have used (|4.31| ) to eliminate the potentials at intermediate stages. 

Now we can combine ( ^4.26 )and ( f4.38| ) to put ( f4.34| ) into a more symmetrical form: 

where 

~ If 1 \ ^ n , muz , , mz . , 

A = ~ 2 (" + u) T8r * + — (b+ + *"> + *T (d - + d+) 

B ± = ±—(l + zd z - Td T ) b± ± — (1 + zdj - Td T ) d± (4.40) 

171 ( 2 t<6 ze ±(j, \ m 

- — z« e +v ^— - —(z — z) cosh0. 

l \ u z J 2 

This is another set of equations in the Lax form, and like ( |4.27|) they behave simply under 
the z «-> ~z, u <-> 6± <->• d± transformations. We can examine the consequences of 
the compatibility conditions dxd z = d z dx and d^d- = d^dy on ( f4.40|) and ( |4.27|) ; this 
generates equations for the potentials which turn to be already implied by the equations 
in Section 4.4.1 and flEggp . 
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4-5. Recapitulation 

Our derivation of the differential equations has been rather circuitous, so it seems 
useful here to present a self-contained review of the final results. The discussion here is 
also designed to be accessible to readers who have skipped Sections 4.1-4.4. 

The main quantities of interest are the correlators r± defined in (|2.20|) . We now 
introduce the variables b±, d± and (f) by 

b± = T—d z hir ± 
m 

d± = ±-bh\nr^ ( 4 - 41 ) 
m 

<j> = ln(r_/r+). 



These variables were introduced earlier as "potentials" , but here we shall simply consider 
( |4.41|) as their defining relations. 

A complete, mutually-independent, set of non-linear partial differential equations 
obeyed by these quantities is 

2 

d x dghiT± = T— e^sinh^ (4.42a) 
d z b + -d z b- = ^(b + + b.) 2 (4.426) 
c\d+ -c\d- = -^(d+ + d-) 2 (4.42c) 
e T<t> (1 + zd z - Td T ) b± - e ±4> (1 + zc\ - Td T ) d± = m(z - z) sinh 0. (4A2d) 



The equations Q4.426D , ( |4.42d ) , and the difference of the two equations in (^.42q|) precisely 



exhaust the set of equations obtained earlier in [0,0. Note that these equations do not 
involve T derivatives, and we have one more equation of this type (one of the two in 
( |4.42a| ) ) than those obtained earlier. The equations in ( |4.42^ ) are new. 

The equations in ( |4.42a| ) and ( 4.42 dp can also be written in the Lax formp3[ 



i.e. as 



the integrability conditions of linear equations obeyed by auxiliary functions dependent 
upon an additional spectral parameter. We need 2 auxilliary functions h±(u), where u 
is the spectral parameter. The linear equations satisfied by h±(u) are ( 4.27|) and ( |4.39| ), 



where we note that on the left-hand-sides of these equations one can replace d 2 — > d z 
d- — > c\ and dy — > (u — 1/u)T8t + (u + l/u)ud u . 
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4-6. Massless Limit 

In this section we will examine the behavior of the non-linear partial differential equa- 
tions in the limit of small m/T. We recall that all correlators are smooth functions of m 
at m = for finite T [21|, ||15| , so the r± obey the following series expansions 



n ~° (4.43) 

r) 1/4 r_=f;(f) 2 " +1 r 2 „ +1 (,T,^). 

n=0 

The prefactor of (m/T) 1 / 4 follows from the overall normalization condition for r± implicit 
in ( |2.22| ), and the behavior of the correlators at m = 0. The coefficients r n in ( |4.43|) are 
dimensionless functions only of zT and zT. Indeed, the r n are clearly spacetime integrals 
over correlators in the massless conformal field theory of two a or \x operators and n 
thermal (energy) operators. These correlators possess a holomorphic/anti-holomorphic 
factorization property, but it is not expected that the result of a spacetime integral over 
them in a cylinder geometry will be in a factorized form. Let us stress here that conformal 
perturbation theory expansions such as (|4.43|) as possible only because T ^ 0, so there are 



no infrared divergences [25 



An interesting question is the convergence of expansions such as ( |4.43|) . Of course, the 
absence of any phase transitions at finite temperature in one-dimensional physical systems 
with short-range interactions implies that there can be no singularity at any finite real 
value of m/T. The only singularities on the real axis consist of powers of exp(— m/T) as 



m/T — > oo (this is clear from the results of [15|). However, phase transitions are possible 
in one- dimensional systems with unphysical complex couplings. In the present case, the 
largest eigenvalue in the spin even sector crosses the largest eigenvalue in the spin odd 
sector when m = inT, and presumably this determines the radius of convergence of Q4.43|) . 

We can now insert (|4.43| ) into the definitions ( [I.41|) and the differential equations 
(4.42.) All the equations organize themselves neatly in series of integer powers of m, and 
demanding that they are obeyed at each order leads to a hierarchy of differential equations 
obeyed by the r n . The zeroth order equations in this hierarchy involve only tq and they 
are 

d z djliiT = 

(l + zd z -Td T )d z lnr = (4.44) 
(l + zdz-Td T )dzlnro = 0. 
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As To is proportional to the spin correlators of the massless, conformal, theory, it is expected 
to factorize into holomorphic/anti-holomorphic parts: 



t (w,w) = V (w)V (w), (4.45) 

where we have introduced w = zT, w = zT. Substituting ( |4.45| ) for tq, we find that 
( [4.44]) are obeyed for any function ^q(w). The form of follows from a, now standard, 
argument based on conformal invariance and its overall scale was determined in |T5 



*oN = 7 • (0 m % , 4.46 

where Ac is Glaisher's constant (In Aq = 1/12— £'(— 1)). Notice that, apart from the scale, 
\&o is determined by the requirement that its only singularities be ^o(w — > 0) ~ m; -1 / 8 
and its periodic images under w — > w + 1/2, and that it have no zeros. 



5. Conclusion 

In conclusion, we would like to stress that the differential equations involving only 
spatial coordinates (the equations (|4.42a-c|) ) would hold for any filling function / in ( |2.22|) . 



This is because these equations follow only from the monodromy relations ( fl.9| ),( f02| ), that 
are insensitive to /. The last equation in (4.42 )would hold for any filling function / that 
depends only on (u + l/u)/T and satisfies /(oo) = 1, /(— oo) = 0. This is because the 
differential equation follows from ( f4.15| ) that in turns follows from ( [1.14| ), together with 



( |4.30|) that depends on the T = behavior of /. This does not mean that the differential 
equations are too general to be useful: they can be quite constraining after imposing 
some limiting behavior of the correlators (like (|4.46[) ) determined by other means (similar 
features are known for the Bose gas 

While we have not derived the differential equations for the boundary problem, we 
can straightforwardly obtain some of them. Indeed, the kernel ( |3.11| ) differs from the zero 
temperature one in [[HJ by the filling function / only. Thus, as for the bulk two point 



function, the differential equations derived in |27] involving derivatives with respect to 
k and x still apply for the non-zero temperature correlator. 

The presence of a non-zero temperature also introduces a number of new physical 
phenomena, not found at T = 0, which are of considerable theoretical interest. Foremost 
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among these is the appearance of irreversible, thermal relaxational behavior in the long- 
time spin correlations (for field-theorists, it is perhaps useful to note here that we are 
referring to correlations in real time; this relaxational behavior has no simple character- 
ization in Euclidean time where correlations are periodic with period 1/T). While this 
relaxational behavior is expected to occur at all non-zero T, there are some important 
distinctions in its physical nature between high T (T ^> m) and low T (T<m), where m 
is mass gap above the ground state of the continuum Ising model. 

(i) T ^> m: The presence of thermal relaxational behavior has been established in this 
limit, and a simple argument based on conformal invariance also allowed exact computation 



of a relaxation rate constant |2q| . In this high T regime, the fluctuations are dominated 
by states with energy Tiuj ~ /c^T, and the spin dynamics are therefore not amenable to a 
description by some effective classical model of dissipative dynamics [[TJJ . 
{%%) T <C m: The physical situation is quite different at low T just above the ordered 
state (we will not comment here on the behavior above the quantum paramagnetic state) . 



Now the most important excitations have an energy Tluj <C ksT [15], and their de-Broglie 



wavelength is therefore much smaller than their mean spacing. This reasoning led to the 



conjecture [jig] that the long-time relaxational spin dynamics is described by an effective 
classical model. 

The study in this paper is a step towards understanding the finite temperature dynam- 
ics of the continuum quantum Ising model, as it crosses over from the high T limit (where 
the dynamic correlations are known in closed form) to its conjectured classical relaxational 
behavior at low T. In all our computations, time evolution is specified by the usual unitary 
Heisenberg operator e lHt / h : where H is the Hamiltonian on the quantum Ising model, and 
there is no coupling to a heat bath. The temperature appears solely in specifying an initial 
density matrix e~ H / T ', which is also stationary under the time evolution. We have taken 
several significant steps towards obtaining dynamic correlation functions under the unitary 
Heisenberg time evolution, and its complete understanding now appears within reach. 
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Appendix A. Direct Approach 



In this appendix, we would like to give a different derivation of the result (|2.13|) . 
To start, we consider the one point function of some operator O: 

In a multiparticle basis for the states \if) >, one has the following resolution of the 
identity when L = oo: 

1 = E h. II dJ W^ IA. • ■ ■ . A.XA.. • ■ ■ . Al- (A-2) 



n=0 



Thus, 



(°)t = \Y,-, fd(3 1 ---d/3 n lf[e- e ^/ T )(Pn,...,Pi\0\(3 1 ,...,f3 n ). (A.3) 

n=0 n ' J \i=l J 

The matrix elements of O are the so called form-factors. As we will explain below, 
one finds that the sum ( |A.3j ) has two kinds of singularities that must be dealt with before 
one can obtain a meaningful expression. The first kind of singularity is a result of working 
with multiparticle states that are defined in infinite volume, i.e. for L = oo. These infinite 
volume singularities manifest themselves as 5(0), i.e. 8((3—(3') with f3 = (3' '. As expected on 
general grounds, these infinite volume singularities are removed upon dividing by Z. The 
second kind of singularity arises from poles in the form factors. These 'pole singularities' 
can be removed by a renormalization of the one-point functions. 

One first needs some general properties of the form factors. We specialize to a theory 
with S-matrix equal to —1. In terms of the operators defined in ( |2.1| ), the multiparticle 
states have the following representation: 

\(3 1 ,...,(3 n ) = Al((3 1 )--.A^(3 n )\0) 

(A.4) 

(#»,..., #l| = (0\A(f3 n )---A\f3i). 

Using the algebra of the operators A, A\ one can evaluate the inner products (f3 n , ...\f3\....) . 

Let A, B denote some ordered sets of rapidities, A = {(3 n , B = {f3[, P' m }, 
and consider the form factor (A\0(x)\B) . If A and B have no overlap, i.e. (A\B) = 0, then 
the form factor obeys the crossing relation (see also next appendix): 

(A\0(x)\B) = (0\O(x)\B,A-i7r), (A.5) 
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where A — in denotes all rapidities shifted by — in. If (A\B) 7^ 0, the form factor is 
expressed as a sum over all ways of breaking up A, B into two sets @ : 

(A\0(x)\B)= (At\0{x)\B 1 )(A 2 \B 2 ) S a , Ai Sb,b,- (A.6) 

A=A 1 uA 2 ;B=B 1 t)B 2 

In this formula, A + denotes rapidities shifted by a infmitesimally small imaginary part irj, 
so that the crossing relation is valid: 

(A + \0\B) = (A + ir]\0\B) = (0\O\B,A-i7r + irj). (A.7) 

The S-matrix factors are defined as the product of S-matrix elements required to bring | A) 
into the order | .4.2,^1), and similarly for B. Namely, 

(A\ = S A , Al (A2,M \B) = Sb,b 1 \B 1 ,B 2 ). (A.8) 

We will refer to the single term in ( |A.6|) with Ai = A, B\ = B as the connected piece and 
the rest as the disconnected pieces of the form factor. 

We will also need the partition function Z. In infinite volume, Z is singular, and it 
is well known how to regulate it on a cylinder of finite length. It suffices for our purposes 
however to work with an unregulated expression for Z, since, as we will see, dividing by 
Z simply cancels similar infinite volume divergences in the sum ( |A.3| ). Therefore, we take 

^ = t^J dJ ^ L (f[»-* /r )ft «' (AJ» 

Up to two particles, using the Faddeev-Zamolodchikov algebra one finds: 

Z = l+ [ d(3e- eiP) / T < p\/3 >+\ ! d/3id/? 2 e- (e(/3l)+e(/32))/T < P 2 ,Wi,fc > + 



2 " (A.10) 
One can explicitly evaluate the lowest order terms in the sum ( |A.3|) , using ( |A.6|) and 
( |A.10|) . One finds that the effect of the disconnected pieces of the form factors is two- fold: 
some of the disconnected terms lead to 6(0) singularities as in Z(T), whereas others lead 
to a modification of the integration measure J d(3. Doing the explicit computation up to 
3 particles, one finds, 

£ e ^ /T w= z((W) 



+ J d[3e-^/ T (l - e~ e ^/ T + e- 2 <^ T ) ((3 + \0\f3) 

+ ^ J dld2 e - {ei+e2)/T (l - e~ ei/T - e- ea/T ) (2+ 1 + |C|1, 2) 

+ dld2d3 e- (ei+e2+e3)/T (3 + ,2 + ,l+|0|l,2,3) + ...) 



(A.ll) 



27 



( dl = dfli, ei = e(Pi),etc.) In the above formula, since we did the computation only up 
to 3 particles, we only verified explicitly the appearance of Z up to the appropriate order, 
depending on which term Z multiplies. 

The above computation, along with some combinatoric checks at higher order, are 
sufficient to understand that to all orders one will find: 

(°>^ = E^ / d(3l (2 n)^ n f[m)e-^ ) ' T {^,...^t\0\^...^n), (A.12) 

n=0 * ^ ^ ' i=l 

where 

= tz — rw ( A - 13 ) 

1 + exp(-e(/3)/T) 

In summary, the non-trivial effect of the disconnected pieces of the form factors is to modify 
the J d(3 integrations by the factor f((3). In the thermodynamic approach, / = p h / (p+p h ). 

The expression ( |A.12|) has another kind of singularity arising from poles in the con- 
nected piece of the form factor. In fact, one of the axioms of the form factor bootstrap 
expresses the residue of the pole in (0\O\...f3, (3' — in, ....) as (3 — > (3' in terms of form factors 
with lower numbers of particles. We will refer to these as pole singularities. In the Ising 
model with S = — 1, for the spin/disorder fields the residue axiom reads 

(0\O\...(3,f3-iir + ir)...) = -(0\O\ ...... .}. (A.14) 

V 



Thus, 



This leads to 



(P+,...,P+\0\Pi,...,Pn}=(^j (0\O\0). (A.15) 



(0) T = hm exp ( [ ^-f(p) e -<^ T ) (O) T=0 . (A.16) 

Note that here we left r\ inside the integration sign, as it may well have to depend on (3 
when one tries to give a more precise meaning to ( |A.16| )(see the next appendix). Here, 
our attitude will be to ignore this (infinite) multiplicative renormalization, our goal being 
to obtain expressions for the two point functions up to an overall normalization, which we 
set somewhat arbitrarily. 

The above procedure leads to meaningful expressions for the two-point correlation 
functions. One begins with 

(O( X )O(0)) T = iV e - £ ,/r»» (A 17) 
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By explicitly evaluating the above sum for low numbers of particles, one reaches the follow- 
ing conclusions. Again, the disconnected terms are of two types, one diverging in infinite 
volume but being cancelled by 1/Z, the other leading to modifications of the integration 
measures as for the 1-point function. We find that the J df3 integrations for the sum over 
states ip and ip' are both modified. To illustrate this, consider for instance the term: 



i J d^^e- e ^/ T ^|O(x)|/3i,^)(^,^|C»(0)|^). 
One of the disconnected pieces is 



(A.18) 



2 



~ I d(3e- E ^/ T (0\O(x)\(3)((3\O(0)\0). 



(A.19) 



This is interpreted as arising from the second term in the expansion of an f(/3) measure 
factor in the sum over states if;'. In addition, the pole singularities can be removed by 
multiplicative renormalization of the fields as described above. Based on this, one finds 
the following expression for the 2-point functions (again, up to an overall normalization): 



WW))* = E is / (um 



|e -e(/9i)/T 



n,m=0 ' \j=l 

m 



X 



H /(#) (0\O(x)\/3[, /3' m , n - in, ft - ztt) 



(A.20) 



\r=l 



x (0\O(0)\(3 1 ,....,(3 n ,f3' rn -in,....,(3[-m}. 

Let us refer to the (3 1 states in ( |A.20| ) as particles and the /3-states as holes. Introducing 
an index e = 1 for particles and e = — 1 for holes, define the particle-hole states as follows: 

\/3i, . . . ,/? n )ei-e„ = \P\ -eiiiT, . . . , /3 n -e n in) 

(A. 21) 

e "- £l (/3 n ,...,/3i| = ( / 9 n + e n z7r,...,/3i + eiZ7r|, 

where e = (e — l)/2. The form- factors of the particle-hole states are the usual ones but 
with the appropriate shifts by in, e.g. 

(O|0(O) |ft , . . . , (3 n ) ei ... en = (0|O(0) |ft - ?iztt, ...,/3 n - e n m). (A.22) 

The expansion ( |A.20|) can the be written as a sum over particles and holes: 

df3\...d(3N 



i r 

(O(*)o(0)> T = Y. m E / 

AT r» ' ^ Li ** 



(27T) 



7V=0 e=±l 

x |(0|O(0)|ft,...,^) ei ... £jv 



i=i 

1 2 



] J (A) exp (-e* {mze pi + mze pi )) 

(A.23) 



where f± are defined in ( |2.12| ). This expression coincides with ( |2.13| ). 
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Appendix B. Derivation of \Bt > 

We work first at T = 0. By considering < 0|O|S >, i.e. in the L channel, and using 



the boundary state of |T9[, we find the expression 



< 0\O(x)\B >= r \n d0 r ■ ■ ■ < 0\O\ - 0!, /3 U . . . , -Pn, P n > 

exp [2ma;(cosh pi + . . . + cosh p n )] . 

(B.l) 

Now, suppose we want to compute this one point function in the crossed, R channel. 
It then reads simply < 0b|(9|0b >, where |0b > is the ground state of the Ising model on 
the half line with a magnetic field at x = 0. 

On the other hand, let us perform crossing on (|B . 1|) . For this, let us first set P = 
i\ — P' . Using the general expression 

< p n+1 ...p n+m \0\p 1 ...p n >=< 0\O\p 1 ,...,p n ,p n+m -iir,...,p n+ i-iir >, (B.2) 
one finds 

< O|0| - px, p U ... , -P n , P n >= (~l) n < -P' n , . . . , -P[\0\p[, ...,p' n > 

(B-3) 

=< — P' n + in, . . . , — P x + i7i\0\P[ — in, . . . , P' n — in > . 

The P' integration runs now from ^ — oo to ^ + oo. Let us move this contour back to 
the real axis - to do so, we request that R has no pole in the physical strip, which will be 
the case for appropriate magnetic field [jT9| . Using that R(P) = R — P), together with 



R(P) = — [R(P — in)]* the one point function ( |B.1| ) can be then rewritten as 

< ~Pn + in, . . ., —Pi + in\C\Pi — in, . . . , P n — in > 
exp {— 2rmo;[sinh Pi + . . . + sinh P n ]} . 



(B.4) 



Expression ( |B.4| ) can be explained by realizing that the boundary changes the nature 
of the Fermi sea of the quantum theory. From ( [B.4] ) , we see that we can represent the new 
ground state \0b > by writing formally as \0b >= ®/3>o {\P — in > —R(P — in) \ — P — in >}. 
Indeed, evaluate now the one point function of O on the half line , < 0s|O|0s >. We first 
have two terms where no R matrix is involved, and which reads naively, using crossing, 

(^ >0 < P\) 0(X) (®f3 >0 \P >) = II Jfl) ' ( R5 ) 
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together with a similar expression with (3 — > — Here is a cut-off necessary to avoid 
the divergences due to particle-antiparticle annihilation. The general form of r\ is easy to 
obtain: at finite L, require that a particle and an antiparticle cannot occupy the same 
state: this leads to a a cut-off in rapidity space such that 5p((3) = mcosh [35(3 = ^, or 
We use this cut-off as the value of rj(f3) in what follows. 

A renormalization of the operator O of course suppresses the term on the right hand 
side of flBjjp . Let us set (®/3>o < (3\) 0(x) (®/3 >0 \(3 >) + (/3 — > -(3) = 1 (much the same 
renormalization would occur in the bulk). Note that there is no x dependence so far. 

Then, there are also crossed terms in < 0s|(9|0s >, which involve the matrix R. 
These are x dependent and cannot be absorbed in a renormalization. Indeed, for a given 
rapidity (3 in the decomposition of the ket |0^ >, in the corresponding bra we can pick 
either (3 or —f3. In that case, compared with (|B.5|) and its (3 — > —f3 equivalent we get a 
relative factor 

<jg> = <-/W> (B 6 > 

<(3+\0\(3> + <-f3+\0\-f3> 2 mLcosh/? ' 

7T 

On the other hand, this particular rapidity where a pairing with a reflected state takes 
place can be chosen with some multiplicity: there are ^ cosh f3df3 states with which one 
can do so in [(3, (3 + d(3] Clearly, the two L-dependent contributions cancel out, leaving 
an overall factor of 1/2. Generalizing to arbitrary number of pairings and using the fact 
that R*((3) = R(—(3), one reproduces the expression (|B.4f ), where the integrals with (3 > 
corresponds to taking the reflected state in the bra, the integral with (3 < to taking the 
reflected state in the ket. 

The foregoing form of \0b > is easily checked (eg by evaluating matrix ele- 
ments between \0b > and any other state) to satisfy the relations |0| A^((3)\0b >= 
R(P)A\-P)\0 B >, using that -[R((3 - ztt)] = ^y. 

To extend the result at finite T, we require similarly for an excited state \(3b > to 
satisfy A{j3)\j3 B >= R*(f3)A(-f3)\f3 B >■ Using R*(f3) = -^y, it follows that an excited 
state is obtained by replacing some of the \f3 — in > —R{(3 — m)\ — (3 + i% > by \(3 > 
+R((3)\ — (3 >. This leads to two types of states in the boundary heat bath \0b,t > 
characterized by a label e = ±1, with relative weights f±((3). As a result, the one point 
function at finite temperature reads then 

1 Here we used the fact that with a boundary, while only positive rapidities are allowed, the 
density is twice as big 
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< -/3 n + e n z7r, . . . , -/?! + eiZ7r|C|/3i - ?iZ7r, . . . , /3 n - e n zyr > 
exp [2mzxei sinh /%] / e . (/?;), 

(B.7) 

where we set e = ^y^. 

One can then get back to the original channel by following the previous transforma- 
tions in the opposite way to find (|3~"3T) , (|3~7l) ■ 

Appendix C. Boundary Field Theory on a Finite Strip 

The techniques of this paper allows us to consider a different geometry, that of an 
infinitely long strip of finite width R. In this geometry one can consider a quantum field 
theory with non-trivial boundary interactions at x = 0, R, where the strip is defined by 
< x < R. Here, R does not have the interpretation as an inverse temperature. Let us 
quantize the theory in a way that views x as the imaginary time. In this quantization the 
Hamiltonian H is the usual bulk hamiltonian with no boundary terms, and the spectrum 
and form factors are the same as for a bulk theory with hamiltonian H. In this picture, 
the boundary interactions are contained in a boundary state as in ( |3.2|) . 

Consider the one-point function of a field 0(x). For x a semi-infinite line with R = oo, 
this correlator was characterized with form factors in [^D[ . Due to the translation invariance 
in the i-direction, this correlator depends only on x. For the finite width strip, one is 
allowed to put different boundary conditions at each end of the strip. Let the boundary 
condition at x = correspond to the boundary state |-B&), and at x = R the boundary 
state \B a ). The partition function in this picture is then 

Z ab (R) = (B a \e- RH \B b ). (C.l) 

The one-point function is defined by 

(B a \0(x)\B b ) R = -\- (B a \e- H ( R -^O(0)e- H *\B b ). (C.2) 
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Since the boundary state involves pairs of particles with opposite rapidity, the states 
which contribute as intermediate states are 



|2n) = | ■,-/?„,/?„). (C.3) 

Inserting two complete sets of states in ( |C.2| ), one on each side of the field O, one 
obtains 



(B a \o iX )\B b ) R = y: — j 1 



(47r) m 



x 



n,m=0 

n \ / m 

yi=l J \i = l 

5' Ol 



(C.4) 



x ((3 n , -p n , . . . , A, -ft | 0(0) | #J, 

where as usual e(/3) = mcosh/3. 

As in appendix A, one finds that the disconnected pieces of the form factors can be 
incorporated by simply modifying the measures for the integrals over (3. Also as before the 
residue singularities lead to a multiplicative renormalization of the fields, which we again 
ignore. It is not hard to see that one obtains 



dp!---d(3 n f°° d(3[---d(3' m 



B a \0( X )\B b ) R = ^ J ——J 

n,m=U 



(4 7r )m 



X 



(C.5) 



vi=l / \i=l 

x ((3+, -/?+ . . . , ft, -[3+\ O(0) | - 0[, ft, . . . , p m ), 

where now 

f(3) = L ^ . (C.6) 

l + e-™<VR b ((3)R*M 

Result ( |C.5| ) can be put in a thermodynamic form by defining f+((3) = f((3) as in 

( |(J.6D , and /-(/?) = 1 — /(/?)• Then, using that R((3 — in) = — J—, one finds 

R((3) 



< B„\0(x)\B b > R = Y. £ jvf / II 



|?r e- 2a:eie(/3i) [^]e(A)/e(A 



(C.7) 



< 0|O| - p 1 ,(3i,... 1 -p N ,(3 N > eit .„, 



This can easily be recovered in a thermodynamic approach. The TBA was developed for 
the partition functions Z a b in |2{J , and one sees that / given in (|C.6|) corresponds precisely 
to the appropriate ratio of densities as determined there. 
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Appendix D. Computations for bulk correlators 



In this appendix we shall provide the missing steps required to establish Q4.30| ) . 
First, we show that the analog of Q4.29|) with T derivatives also holds. As in Q4.16| ) 
we have 

T<9 T hi r± = ±Tr((l-R±)Td T K). (D.l) 



From the defining relations (|2.34| ) and (|2.35| ) we have 



TP) \C ( \ IMM xCf ,Td T f(v) 



2/(«) 



2/(«) 



(D.2) 



Inserting this into ( 4.16 ), and using Q4.1| ) we get 



Td T lnr± = Tr R±(u,u 



Td T f(u) 
/(«) 



(D.3) 



where the symbol Tr implies an integral from — oo to oo over variables u 7 v, . . . that appear 
in its argument. Now using the explicit expression for R in ( f4.4|) and ( |4.7|) we get 



t« i fa f \ f \ ( \a ( \ jl 9-(u)g+(u)\ Td T f(u) 
Td T lnr ± = -Tr nd u g-{u)g + (u) - g-{u)d u g+{u) ± 1 



(DA) 



In particular, we have 



T9 T In r_/r + = -Tr . 



(D.5) 



In the following, we will show that the right hand side of ( p.5| ) is simply related to the 
potentials which were introduced earlier. We have also examined T8t ln(r_r_|_), but its 
value does not appear to be expressible solely in terms of the potentials of Section 4.3. 

Apply the operator (1 — R±)T<9x to both sides of the first equation in ( |4.6| ). Using 
( p.2| ) and simplifying we get 



Tdrfju) 



dvTl±{u, v ) — g±{y) . 



(D.6) 



Also, from (ET34D we have 



TP, ( \ Td rf( u ) i ^ 

Tdre{u) = ~vW e{u) - 



(D.7) 
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We are now ready to take T derivatives of the potentials. We have from the definition 
( CT) and (ED, flD3) 

Td T c± = Td T Tr(ed_i#±) 

= IV (ed_i(l - R±)^y4±) (D g) 

= (l T c ± )Tr^ T d_ 1 ^^ ± 



where in the last step we have used the intertwining relations ( |4.10| ). Finally, using the 
parametrization ( f4.21|) , and comparing with ( p.5|) , we get 



Td T (ln(r_/r+) - 0) = 0. (D.9) 

The combination of ( f4.29| ) and ( p.8| ) implies that t_/t_|_ = Ce^ with C a pure number 
independent of z, z, and T. We now determine C by computing the large z, z behavior of 
t± and <f) at T = directly from the form factor expansion. Indeed, simply by evaluating 
the first two terms in the summation in (|2.22j ) we have 

t± = 1 ± 7 + . . . T = 0, r -> oo (D.10) 
\J1ixmr 



where r 2 = 4,2 z = x 2 + t 2 . Inserting (|D.10| ) into the first differential equation in Q4-28 ) we 
get 

2g— mr 

+ ... T = 0, r ^oo (D.ll) 



Finally, inserting ( p.lO| ), ( p.ll| ) into t_/t+ = Ce^ and matching the large r behavior, we 
get C = 1. 
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